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We consider heavy quark production in high energy pA collisions and investigate the 
contribution of interactions of valence quarks of proton with the nucleus. The often made 
assumption that valence quarks of proton can be factored out is justified only if the nucleus 
saturation momentum is much smaller than the heavy quark mass. This is not the case 
in phcnomenologically relevant situations. Breakdown of factorization manifests itself in 
substantial decrease of the cross section at large total and small relative transverse momenta 
of the heavy quark - antiquark pair. 

I. INTRODUCTION 

The most general expression for heavy quark pair production in high energy pA collisions was 
derived in [1]. It takes into account interaction of valence quark q v , intermediate virtual gluon and 
the produced qq pair with the heavy nucleus target. That expression is quite bulky and unfriendly 
for numerical calculation because it involves multi-dimensional nested integrals over the oscillating 
integrands. Therefore, one usually considers an approximation in which interaction of valence quark 
with the nucleus is neglected leading to collinear factorization of the gluon distribution function 
on the proton side [2-7]. This approximation turns out to be valid only in the limit of very heavy 
quarks m S> Q s , where Q s is the saturation momentum of the target nucleus. Since, Q s turns 
out to be of the same order of magnitude as the charm and bottom quark masses, interactions of 
valence quarks must be taken into account. It is the goal of this article to develop an approximation 
that takes into account interaction of valence quark with the nucleus and still allows one to make 
relatively simple numerical analysis of heavy quark production. 

To develop a consistent approximation we need to ensure the gauge invariance of each term in 
the expansion. The derivation is presented in the main part of this article, which is structured as 
follows. In Sec. II we review the result of [1]. In Sec. Ill we discuss the factorization limit of these 
general formulas. Contribution of valence quarks is calculated in (IV). Our main result is given 
by Eqs. (78), (67), (75). The effect of valence quark interactions on double inclusive cross section of 
charm production is shown in Fig. 3. We denoted the transverse momenta of quark and antiquark 
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as ki and k2, the fraction of the gluon's light-cone energy carried by the heavy quark as z, the pair 
total momentum as q = k\ + k,2 and their relative momentum as t = (l — z)ki — zk2, which is related 
to the invariant mass of the pair as M 2 = (m 2 + (?) / z(l — z). It is seen on this figure that deviation 
from the collinear factorization of proton is negligible at small q, which corresponds to valence quark 
being collinear with the valence quark, and large £, which corresponds to large invariant mass of the 
pair. However, it grows as q increases and £ decreases. It appears that the collinear factorization of 
proton is a reasonable approximation for total cross sections which receive the main contribution 
from low transverse momenta. However, for single and double-inclusive spectra it overestimates 
the cross section by orders of magnitude. 

II. HEAVY QUARK AND ANTIQUARK PRODUCTION IN PA COLLISIONS: 

GENERAL RESULT 
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(1) (2) (3) 

FIG. 1: Diagrams contributing to quark-antiquark pair production in the light-cone gauge. (1) Incoming 
valence quark emits a gluon, which splits into a quark-antiquark pair before the system hits the target. (2) 
Valence quark first emits a gluon, after which the system rescatters on the target nucleus, and later the 
gluon splits into a quark-antiquark pair. (3) Valence quark rescatters on the target nucleus, after which it 
produces a gluon, which splits into a quark-antiquark pair. 

The diagrams contributing to quark-antiquark pair production in the light-cone gauge are shown 
in Fig. 1. The vertical dashed line depicts interaction of the projectile partonic system with the 
target. Since the production time of the partonic system is proportional to the collision energy, it 
is much larger than the interaction time, which is taken to be vanishingly small. The calculation is 
done in the light-cone perturbation theory [8], along the lines outlined in [10]. It is convenient to 
work in coordinate space where the diagram contributions factorize into a convolution of Glauber- 
Mueller multiple rescattering [9] with the "wave function" parts, which include splittings q v —> q v g 
and g —> qq. Since eikonal multiple rescatterings do not change the transverse coordinates of the 



3 



incoming quarks and the gluon, we can calculate the "wave functions" in transverse coordinate 
space by calculating the diagrams in Fig. 1 without interactions. We denote momenta of the out- 
going quark and anti-quark as k\ and k 2 correspondingly. We assume that the gluon is much softer 
than the proton, i.e. k± + + k 2+ <C p+ k 2+ *C p+, where p + is the typical light cone momentum 
of the valence quarks in the proton. Expressions for the wave functions in the momentum space 
can be found in [1]. The light cone "wave- functions" in transverse coordinate space are defined as 

<U*i,* 2 ;.) = J^^e-^^V^ikut*), , = 1,2,3, (1) 

where the superscript (i) corresponds to one of three diagrams in Fig. 1. Here we assume that the 
transverse coordinate of the valence quark, which emits the gluon (which splits into a qq pair) is 
and denoted the transverse coordinate of quark and antiquark by X\ and x 2 correspondingly. 
The quark-antiquark production cross section is proportional to the product of the sum of the 
light-cone "wave-functions" in the amplitude and sum of the light-cone "wave-functions" in the 
complex-conjugated (c.c.) one, averaged over the quantum numbers of the initial valence quark 
and summing over the quantum numbers of the final state quarks. Since transverse momenta of 
quark k\ and antiquark k 2 are fixed, their coordinates in the amplitude and in the c.c. amplitude 
are different. We will denote the corresponding coordinates in the c.c. amplitude by x[ and x' 2 . 
The resulting quantity appearing in the cross section is given by 

<f> ij (x 1 ,x 2 ;x' 1 ,x' 2 ;z) = — ^ a ,(x 1 ,x 2 ; z)^*,(x' 1 ,x' 2 ; z) , t,j = 1,2,3. (2) 

C cr,cr',a,b 

Here the sum over gluons' colors a and b simply implies a calculation of the color factors of the 
relevant diagrams, including traces over fermion loops. 

The double-inclusive quark-anti-quark production cross section in pA collisions in the quasi- 
classical approximation, i.e. neglecting quantum evolution of the partonic system, 



d 2 kid 2 k2 dy dz d 2 b 



= -^— R [ d\ 1 d 2 x 2 d 2 x[d\' 2 e- tk ^~ x '^ k ^~ x '^ 
4(2^)° J 

3 

x X] ^'(^I'^s;^!,^;^) Sij (xi,x2 ; x[,x 2 ;z) . (3) 
»J=1 



Here y is the rapidity of the s-channel gluon, which splits into the qq pair. are the 5-matrix ele- 
ments corresponding to the i'th wave function in the amplitude and j'th in the complex-conjugated 
amplitude. Since the quark and the anti-quark are most likely to be produced close to each other 
in rapidity, one can think of y as the rapidity of the quarks, b is the impact parameter of the 
proton with respect to the nucleus. 
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The single inclusive quark production cross section is obtained from (3) by integrating over one 
of the quark's momenta and multiplying by 2 to account for both quarks and anti-quarks: 

da q ^ qX 



d 2 kdyd 2 b 2(2 



^—j [ d 2 x 1 d 2 x 2 d 2 x' 1 C dze- ik ^- x 'J 
'tt) J Jo 



x ^2 $ij (xi,x 2 ;x' 1 ,x 2 ;z) Si j (x 1 ,x 2 ;x' 1 ,x 2 ;z) , 



(4) 



where y is the rapidity of the produced (anti-)quark. 

Explicit expressions of the "wave-functions" can be found in [1]. Upon summation over A they 
read 



2g 2 T a T b 



F 2 (x l7 x 2 ; z) — [(1 - 2z)x 12 ■ u + i a m x l2j ] 5 a y 

X\ 2 u 



^L(x 1 ,x 2 ;z] 



tf^Uasi, ^2; z) = -^(xi, a; 2 ; z) - *^,(*i, * 2 ; z) , 



+ Fi(xi,x 2 ; z) -am(u x + iau y ) <5 CTi 

- 2 5 CT)CT / z(l - z)F (x 1 ,x 2 ;z) 
2g 2 T a T h 



(2nf 



mKi(mxi 2 ) — — ^ [(1 - 2z)xi 2 ■ u + i a UiX 12j ] 5 a y 



+ K (mx 12 ) -^am(u x + i a u y ) b c - a i 



u 



(2) 



(5) 

(6) 
(7) 



where €±2 = 1 = —621, en = 622 = 0, and, assuming summation over repeating indices, CijUiVj = 
Also xi 2 j denotes the jth component of the vector x\ 2 . Transverse coordinates of 
gluon in the amplitude is 



U = ZX\ + (1 — z) x 2 . 



(8) 



We denote u = x\ 2 = x\ — x 2 , x\ 2 = \x\ 2 \. Similarly, gluon's transverse coordinate in the c.c. 
amplitude is u' = zx\ + (1 — z)x' 2 with u' = \u'\ and x' 12 = x\ — x' 2 , x' 12 = \x'i 2 \. To perform the 
Fourier transform of (1) we first introduce the following auxiliary functions 

F 2 (xi,x 2 ;z) = J dq J\ (qu) K\ [x\ 2 \J m 2 + q 2 z(l - z)^j \J m 2 + q 2 z(l - z) , (9) 

F 1 (x 1 ,x 2 ;z) = J dq Ji (qu) K (x 12 m 2 + q 2 a(l - zf\ , (10) 

F (xi,x 2 ;z) = J dqqJ Q (qu) K (x 12 ^ m 2 + q 2 z(l - z)^j , (11) 
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where q = k\ + fc 2 . Substituting into (2) we derive 

<f> 11 (xi,x 2 ;x' 1 ,x 2 ;z) = 4C F (— { F 2 (x 1 ,x 2 ; z) F 2 (x' 1 ,x' 2 ; z) \ 4(1 - 2zf 

V vr / I x\2X l2 uu' 



x (x 12 ■ u)(x' 12 ■ u) + (eijUiX 12 j) (ekiu' k x' 12 i)] + F 1 (x 1 ,x 2 ; z) F 1 (x' 1 ,x 2 ; z)m 2 
+4z 2 (l - zf F (aJi, x 2 ; z)F (x' 1 , x' 2 ; z) - 2z(l -z)(l- 2z) 



u • u' 



xyi-u 

F 2 {xi,x 2 ; z) 

X\ 2 U 



X 12 ' U 



x F (x[,x 2 ;z) + 7 , F 2 (x[,x' 2 ;z)F (x 1 ,x 2 ;z) 



x 12 w 



(12) 



•$> 22 (x^x 2 ;x\,x' 2 ;z) = AC F (—\ m 2 \ Kiimx^K^mx'u) / 2 /2 [(1 - 2z) 2 

'It * ?x' 1 

x (xi2 • m)(»'i2 • «') + (tijUiXi2j){eklu' k x' 12l )] + K (mx 12 )K (mx' 12 )^^ \ > ( 13 ) 

$i 2 (xi,x 2 ;xi,x 2 ;2:) = -4Cf( — J m|F 2 (xi,x 2 ;2;) J FCi(mxi 2 ) r ^[(1 - 2z) 2 

V TT y [ zi 2 ar 12 uu' 2 

x (xi2 • u) (x' 12 ■ u') + (eijUiX 12j )(e k i u' k x' 12l )] + mF 1 (x 1 ,x 2 ; z)K (mx' 12 ) ^ 2 



x' • u' ~) 

-2z(\ - z)(l - 2z)^^ 2 -F (x 1 ,x 2 ;z)K 1 (mx' 12 ) \ , (14) 
x 12 u J 

All other products of the wave functions can be found using relations 

®33(xi,x 2 ; xi, x' 2 ; z) = $n(xi, x 2 ; xi, x' 2 ; z) + $ 22 (xi, x 2 ; xi, x' 2 ; z) + $i 2 (xi, x 2 ; x[, x 2 ; z) 

+ $ 2 i(xi,x 2 ;x' 1 ,x 2 ;z) (15) 
$i3(£Ci,x 2 ;xi,x 2 ;2;) = - $n(xi, x 2 ; xi, x 2 ; 2) - $i 2 (xi, x 2 ; xi, x' 2 ; z) (16) 
$ 2 3(a:i,x 2 ;xi,x 2 ;2:) = - $ 2 i(xi, x 2 ; xi, x' 2 ; z) - $ 22 (xi, x 2 ; xi, x' 2 ; z) (17) 

and 

$y(xi,x 2 ;xi,x' 2 ;z) = <Z>* i (x' 1 ,x' 2 ;x 1 ,x 2 ; z). (18) 

Here Eqs. (15), (16), (17) follow from (7). Eq. (18) allows one to obtain $21, $31 and $ 32 from (14), 
(16) and (17). 

Rescattering of q v , q v g and q v qq configurations on a large nucleus brings in different factors, 
which we label by Sij. For the case of single-quark inclusive production cross section they were 
calculated in [3], while in the case of double-inclusive production in [7, 11]. These factors can 
be written down as a combination of color dipole and color quadrupole scattering amplitudes. 
For notational simplicity we will assume the large-iV c limit in which only color dipoles survive. 
Generalization of our results beyond the large N c limit is straightforward, though very bulky. 



Sij's are conventionally expressed in terms of the gluon saturation momentum defined as 



Q 2 s =AW sP T{b) 



(19) 



where p is the nucleon number density in the nucleus and T{b) the nuclear profile function. The 
scattering factors (which are proportional to the two-point correlations functions) are given by 



Sf(x) = e 



-±x 2 \n{l/xk)Q 2 a 



S A {x) = e 4 



±x 2 \n(l/xk)Q 2 3 



for quark and gluon color dipoles correspondingly. Using these definitions we have [1] 



S 1 i(xi,x 2 )x' 1 ,x' 2 ;z 
S22(xi,X2]x' 1 ,x' 2 ;z 
S33(xi,X2;x' 1 ,x' 2 ;z 
S 12 (x 1 ,x 2 ;x[,x' 2 ;z 
S 2 3(x 1 ,x 2 ;x , 1 ,x 2 ;z 
Si 3 (x 1 ,X2)x' 1 ,x 2 ;z 



Sf(xi - x[)Sf(x 2 - x 2 ) , 

Sa(u - u') , 

1, 

Sf{xi - u')Sf{x 2 - u') , 
S A {u) , 

Spix^Spix^ , 



(20) 



(21) 
(22) 
(23) 
(24) 
(25) 
(26) 



where A is an infrared cutoff. All other S^'s can be found from the components listed in (21) using 



Sij(xi, x 2 ; Xi, x 2 ; z) — Sji(xi, x 2 ; x%, x 2 ; z) 



(27) 



similar to (18). 



III. LIMIT OF PROTON COLLINEAR FACTORIZATION 



A. Gluon scattering on heavy nucleus 



q, u 





O 



FIG. 2: Gluon-nucleus interactions contribution in the proton collinear approximation. 



Suppose now that valence quarks do not contribute to the quark-antiquark production. In 
this case the relevant diagrams in the light-cone perturbation theory are depicted in Fig. 2. The 
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corresponding "wave-functions" in the two cases differ only by the relative sign and are given by 

KZtoM- .-J? g^W^g>. (28) 



Since integration over q yields (2n) 2 6(u) we can consider the light-cone wave functions in the 
reference frame where q = and integrate only over the remaining momentum q' = k,2 = —k\. 

x |e A • [(1 - z)k! - zk 2 ](l -2z + Xa)5 a y + ^=<rm(l - A<t)<5 ct ,_*'} (29) 
It is convenient to change variables from x± i2 and x' 12 to u, u', x\2 and x' 12 - 

Xi = u + (1 — z) Xyi x'i = u + (1 — z) x' 12 (30) 

X2 = U — Z X\2 x' 2 = 11 — Z x' 12 ■ (31) 

Performing Fourier transformation yields 



^ 9 {l q l{xi,X2) = ^—5(u)\imKi(mxi2)- — — — (1 - 2a + \a)8 C(T > 

-5 a - a -^=am(l - Xa)K (mxi2)\ . 



(32) 



Multiplying by the contribution of the c.c. diagram, summing over the quantum numbers of the 
final states and averaging over the quantum numbers of initial particles according to 

^(x 12 ,x' 1 2,z) = ^^E^V?(*i.^S5?(*i'.«2) (33) 

we derive 

^(x 1 2,x' 1 2,z) = ^5(u)5{u')m 2 { Xl2 ' X ^ [(l - zf + z^K^mx^K^mx'^) 
[Air) i xi2X 12 

+K (mxi 2 )K Q (mx' 12 ) j . (34) 

In particular, in the chiral limit (34) reduces to 

4>^(xi2,x' 12 ,z) = ^SMSW^ll^PM , (35) 

\ z7r ) X \2 X Y2 

where 

P qg (z) = \[{l-zf + z 2 ) (36) 

is the the splitting function. 
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The cross section reads 
fog-yqqX 



d 2 hd 2 k 2 dydz (2 



x^9(x 12 ,x' 12 ,z) 



K-Mxp, q 2 ) J" I d 2 x l d 2 x 2 d 2 x' l d 2 x l 2 d 2 u d 2 v! e -i*i-(*i-*i)-iM*2-^) 



- x\)S F {x 2 - x' 2 ) - Spix'JSFix't) - S F ( Xl )S F (x 2 ) + 1 



, (37) 



where S± is the transverse cross-sectional area of nucleus. Integrals over u and u' are trivial due 
to the delta-functions in (35). It is convenient to express the final result in terms of the transverse 
moments q and I defined as 

q = k 1 + k 2 , £ = (1- z)ki - zk 2 . (38) 

q is the gluon transverse momentum, i.e. total momentum of the quark and antiquark. I is the 
relative transverse momentum of the pair. Invariant mass of the pair M can be written in terms 
of I as follows 

M 2 = (h + k 2 ) 2 = . (39) 

Substituting (30) and (31) and performing integration over the impact parameter b' = (x^ +x' 2 )/2, 
which yields S±, we get 

cPicPqdydztfb "(2tt) 4 tt^'' 7 > J d Xl2Ci x ^ e 

x {K 1 (mx 12 )K 1 (mx' 12 ) X ^^[(l - zf + z 2 } + K {mx 12 )K {mx' 12 ) 1 



(40) 



S F ({1 - z)(x 12 - x' 12 ))S F (z(x 12 - x' 12 )) 
- S F ((1 - z)x' 12 )S F (zx' 12 ) - S F ((1 - z)x 12 )S F (zx l2 ) + l| . 
This is the formula derived before in [2, 3]. 

B. Approximation of the general result of Sec. II 

Now we would like to find an approximation to the general formulas of Sec. II that lead to 
the same result (40). Collinear factorization of proton involves several assumptions. First, the 
characteristic transverse momentum scale of proton A is assumed to be much smaller than the 
produced quark mass. As a consequence, the size of color dipoles u and u' is much larger than 
that of X\ 2 and x' 12 i.e. 

212 < u < A -1 . (41) 



Second, factorization requires that gluon coordinate be the same in the amplitude and in the c.c. 
one independently of the produced quark and antiquark coordinates. This is ensured if 



\u -u'\<^Q s l 



(42) 



Taking approximation (41) and (42) in (21) leads to the following expressions of the scattering 
matrix elements: 



S 11 (x 1 ,x 2 ;x' 1 ,x' 2 ;z 
S22(x 1 ,x 2 ;x' 1 ,x' 2 ; z 
Sss(xi, x 2 ; x\,x' 2 \ z 
S 12 (x 1 ,x 2 ;x' 1 ,x 2 ;z 
S 2 3(xi,x 2 ;x' 1 ,X2;z 
Si3(xi,X2;x' 1 ,x' 2 ;z 



S F {{1 - z)(x u - x' 12 ))S F (z(x 12 - x' 12 )) 

1, 
1, 

S F ((1 - z)x 12 )S F (zx 12 ) , 
S A (u) , 
S A (u) . 



(43) 
(44) 
(45) 
(46) 
(47) 
(48) 



Double inclusive cross section (3) can be written as 

d(J q v ^qqX 



(PhtPhdydztPb "4(2tt)6./ U "" " " •' r "" ■' 12 



d 2 U d 2 u'd 2 X 12 d 2 X ' l2 e -i(kl+k 2 y(u-u') e ^[(l-z) kl -zk 2 H Xl2 - X ' 12 ) 



, x 2 ; Xi, x 2 ; z) Sij(x±, x 2 ; Xi, x 2 ; . 



(49) 



To take integrals over u and u' we note that that Sij, with i,j / 3 do not depend on these 
coordinate use the following identities 

f e- iq - u F 2 ( Xl ,x 2 - z)-d 2 u = -1-Ki\K x (xviyJrr? + zjl - z)q 2 ) v / m 2 + zjl - z)q 2 , (50) 
J u q 

[ e-^F^x^xr, z)-d 2 u = -2vri-|^o( W™ 2 + z(l - z)q 2 ) , (51) 
J u q 

J e- iq - u F ( Xl ,x 2 ; z)d 2 u = 2nK {) {x 12 yJm 2 + z{l - z)q 2 ) , (52) 



/ 



e~ qu ^ d 2 u = -2m 



q 



(53) 



In derivation of these identities we used the orthogonality of the Bessel functions 

/ uJ a (qu)J a (qu)du=-5(q-q). 
JO 1 



(54) 
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Denoting s = \Jm 2 + z(l — z)q 2 we write 
J d 2 u J d 2 u' § ll {x l ,x 2 -x' 1 ,x' 2 -z)e~ iq - u+iq - u ' 

= 4C F ( (2vr) 2 [k^x^K^x'^s 2 - 2z) 2 (x 12 • q)(x' 12 ■ q) 

+ (.eijqixnj) {e k iq k x' l2l )} + Kq{xi 2 s)Kq{x' 12 s)-^- + 4z 2 (l - z) 2 K (xi 2 s)K (x' 12 s) 



- 2z(l - z)(l - 2*) 



(-i) K 1 (x 12 s)K (x' 12 s)s H ; — T i K 1 (x' 12 s)K (x 12 s)s 

xi 2 q z x\ 2 q l 



-iq-u+iq-u 



(55) 



J d 2 u j d 2 u $ 22 (xi,x 2 ;x' l ,x' 2 ;z)e %q 

= 4C F ( — ) (2i:) 2 m 2 {K l {x l2 m)K l {x l l2 m) \—A{\ - 2z) 2 (x l2 ■ q)(x' 12 ■ q) 

+ {(-ijQiXuj) (e k iqkx' 12 i)} + K (£i 2 m)ifo(^2 m )^2 j > ( 56 ) 

y d 2 u y dV^tJi,^;^,^;^^'" 1 ' 

= -4C F ( — ) (27r) 2 m{K 1 (x 12 s)K 1 {x' 12 m)s ^[(1 - 2z) 2 (cc 12 • q)(x' 12 ■ q) 

\ 7T J { x\ 2 x 12 (t 

+ (eijqixwj) (ekiqkx'm)] + mKo(x 12 s)K {x' 12 m)^ 

- 2z(l - z)(l - 2z) 3C ^iK 1 (x' 12 m)K (x 12 s)\ . (57) 

x 12l J 

In the factorization limit, only terms with i,j = 1,2 contribute as we demonstrate in the next 
section. Therefore, in this section we concentrate only i,j = 1,2 terms. The collinear limit comes 
from the logarithmically enhanced terms that correspond to the most singular in 1/q behavior. 
Thus, in order to obtain the collinear factorization of proton, we need to assume that 

q-^2m, (58) 

which implies that \u — u'\ 3> mT 1 . In view of (41) it means that we require Q s <C m. Using (58) 
in Eqs. (55), (56), (57) we conclude that they coincide (apart from the relative sign): 

y d 2 u j d 2 v! $n e -w-«+w-«' ~ j d 2 u j d 2 u' $22 e - iq u+iq u ' « - Jd 2 uJ d 2 u' $i 2 e - iq u+iq u ' 
« 4C F ( — ) (2vr) 2 m 2 { K^x^K^m) [(1 - 2z) 2 (x 12 ■ q)(x' 12 ■ q) 

\ 7T / ^ «12^l 2 Q 

+ (eijQiXuj) (ekiqkx' 12l )} + K (x 12 m)K (x' 12 m)^ . (59) 
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Substituting (59) and (43)-(48) into (49) and averaging over directions of q using (qiqj) = q 2 5ij/2 
and eijeu = 25 ji we obtain 



dH d 2 q dy dzd 2 b~ (2vr)4 ^ ^ J g2 j d x 12 d x 12 e 



J d 2 x\ 2 d 2 x\ 

x ( K 1 (x 12 m)K 1 (x' 12 m) Xl ^^ [(1 - zf + z 2 } + X (ii2m)Jro(V)l 
x |sf((1 - z)(xi2 - x' 12 ))S F (z(x 12 - x' 12 )) + 1 
- S F ((1 - z)x 12 )S F (zx 12 ) - S F ((1 - z)x' 12 )S F (zx' 12 )} . 

Unintegrated gluon distribution function (p(x p , q) denned such that 

rQ 2 

XpG(x p ,Q 2 )= / (p(x p ,q 2 ) dq 2 . (61) 



(60) 



Since in the leading logarithmic approximation 



xG(x,Q 2 ) = a ^ln^, (62) 



we can write (60) as 

^-^m 2 ^x p ,q 2 ) f d 2 x 12 d 2 x' 12 e-^-<^ 



acJ fact 



d 2 £d 2 qdydzd 2 b (2 7 



x \K 1 {x 12 m)K 1 {x\ 2 m) '^^ [(1 - z) 2 + z 2 ] + XoO^m^oO^™) 

x js>((l - z)(a?i2 - x' 12 ))S F (z(x 12 - x' 12 )) + 1 

- S F ((1 - z)x 12 )S F (zx 12 ) - S F ((1 - z)x' l2 )S F (zx' l2 )} . 



(63) 



As expected this coincides with (40). 

If we are interested in single inclusive cross section, we have to integrate (63) over q and multiply- 
by 2, which yields 

&> vi^ z I ^ I ^g(*,i/i* 12 -*' 12 i) .-*<■■"" 

x (^(mxi^iKjj^^Kl - z) 2 + z 2 } + if (ma;i2)iiro(mxi 2 ) 1 
I xi 2 x 12 J 

x |sf((1 - z)(xi2 - x' l2 ))S F (z(xi 2 - x' l2 j) + 1 



- S F ((1 - z)x 12 )S F (zx 12 ) - S F ((1 - z)x' 12 )S F (zx' 12 )} . 



(64) 
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IV. CONTRIBUTION OF VALENCE QUARKS 

In the factorization approximation valence quarks do not contribute. Indeed, this is evident 
upon substitution (59) into (15)-(17): $13 = $23 = $33 = 0. Contribution of valence quarks 
emerges, along with other contributions, when we keep the logarithmically sub- leading terms in 
(55)-(57). In other words, we are still going to work within the approximations (41), (42), but relax 
the approximation (58). Thus, contribution of valence quarks arises when Q s > m. 

It is convenient to rewrite the last line in (3) using (15)-(17) as (dropping the summation sign) 

SijSij = (QijSij); + (QijSijijj (65) 

where 

(QijSij^ =$11-Sll + $l 2 S*12 + $21^21 + $22522 > 
($^)„ =$ll(S 3 3 " ^23 " S 31 ) + $1 2 (<S 33 " 5i 3 " S 32 ) 

+ $2l(5'33 — 5*23 — 531) + $22(5*33 — ^32 — S23) • (66) 

Because of (41), (^ijSij) l depends only on x± 2 and x' 12 , whereas ($jj<%) n depends only on u 
and u' . We will denote the corresponding contributions to the cross section by as o\ and o\\ 
correspondingly. We detailed the calculation of o\ in Sec. IIIB. The double-inclusive cross section 
averaged over the directions of q is 



da qv->qqX 1 f a s \ 2 f 2 ,2/ -ie-(x l2 -x') 

ndqWedydztPb-JW F WJ J U 

x 



ai2 'f 12 [(l - zf + z 2 ]^[K 1 {x 1 2s)K 1 {x' 12 s)s 2 S 11 + K l (x 1 2m)K l (x' l2 m)m 2 S22 

X 12 X \2 Q 

2 

- K 1 (x 12 s)K 1 (x' 12 m)msSi2 - K 1 (xi 2 m)K 1 (x' 12 s)msS 2 i] + [K (x 12 s)K (x' 12 s)S n 
+ K (x 12 m)K (x' 12 m)S 22 - K (x 12 s)K (x' 12 m)Si 2 - K (x 12 m)K (x' 12 s)S 2 i] 



+ 4z 2 (l - z) 2 K (x 1 2s)Ko(x , 12 s)S n ^ . 



(67) 



Integral over q is logarithmically divergent in the UV, with the cutoff Q such that q < Q. 
Now we turn to the contribution an- Substituting (43)-(48) we get 

= d 2 ud 2 u' e' iq - [u - u>) / d 2 x^ d 2 x\, e -^-^-^ 2 ) 

d?k 1( Pk 2 dyd?bdz 4(2vr)6 J auaue J ax 12 ax 12 e 

x [1 - Sa(u) - S A {u')} {^ii(xi,X2;x' 1 ,x , 2;z) + $12 (an, x 2 ; x[, x 2 ; z) 

+ $2i(a5i,£c 2 ; a3i,£c 2 ; z) + $22(0:1, x 2 ; x[, x' 2 ; z)} . (68) 
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Since the scattering matrix elements are independent of x± 2 and x' 12 , we can integrate the wave 
function over these variables. The integrations are performed employing the following identities 

1 



/ 



d 2 x 12 e~ ie - Xl2 F 2 (x 1 ,x 2 ,z)^ = -27ri£ dqJ 1 (qu)—— {j 

xu J I 1 + z(l - z)q z + m z 



2txH 



(£ 2 + m 2 )u 



l 2 + m 2 f £ 2 + m 2 



2ml 



{I 2 + m 2 )u 



h(u,£,z). (69) 



J d 2 x 12 e- ie - xl2 F 1 (x 1 ,x 2 ,z) = 2vr^ dgJi(gu)- 



+ z(l — z)q 2 + m 2 



2vr 



{I 2 + m 2 )u 



fi(u,£,z) . 



(70) 



J d 2 x 12 e i£ ' Xl2 F (x ll x 2l z) =2ir J dq J (qu) 



i 2 + z(l — z)q 2 + m 2 



2tt 



-Kn U\ 



I 2 + m 2 



2vr 



z(l — z) \ V z d — z ) / ^ ~~ z 
Using these formulas we have 



-f (u,e,z) . 



J d 2 x 12 J d 2 x' 12 e- it - x ^ +ii <^ 11 (x 1 ,x 2 ;x' 1 ,x' 2 ;z) = 

= ±C F (^) 2 (^) 2 { J U , 2 (£2+ 1 m2)2 [( 1 " 2z f( u ■ W ■ 4 + 



x fi(u,t,z)fi{u ,t,z) + m 



,u ■ u 



+ 4z 2 (l-z) 



2z(l - z)(l - 2z) 
u' l 1 



z 2 (l-z) 2 
ul 



u 2 u' 2 (I 2 + m 2 ) 2 
fo(u,£,z)fo(u',£,z) 



(-0 



z(l — z) ^ 2 + w. 2 



/i(n,£, z)fi(u',t,z) 



Mu,e,z)Mu',e,z) 



,'2 



(1-z) £ 2 + m 2 



/o(n,£, z)fi(u',£,z) 



(71) 



(72) 



d #12 y d x 12 e lt ' 



+ 



AC F [ — ) (2n) 2 m 2 , , , , 

V 7r / [ m z u z u' z (i z + m 2 ) 

77, ■ (7,' 1 



^22(xi,x 2 ;x' 1 ,x 2 ;z 
1 1 



2 [(1 - 2z) 2 {u ■ i){u' ■ I) + {e^Uitj^u'^)} 



u 2 u' 2 (t 2 + m 2 ) 2 



(73) 
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J d 2 x 12 J d 2 x[ 2 e- u - x " +i£ - x '™$ 12 (x 1 ,x 2 -,x' 1 ,x' 2 -,z) = 

. j 

x fi(u,£,z) + m^ ri 2— 2- ^fi{u,l,z) 

- mi - «)d - 3) P + 1 m , /°("'.^)} p") 

Fourier transforms of other wave functions can be derived using (15)— (17). Substituting in (68), 
the resulting expression averaged over the directions of I (to make it less bulky) reads 



7rd£ 2 d 2 qdyd 2 bdz 4(2vr) 4 r \ vr 

u - u' 1 



X 



;(V , (f , m 2 ) 2^(( 1 " *)* + **) +m 2 ](/ 1 (u,£, 2! ) - l)(/i(u',*,z) - 1) 
+ 4f (u,e,z)f (u',e,x)} [1 - S A (u) - S A (u')] . (75) 

To obtain the single inclusive cross section we can integrate over I or over q and multiple by 2 
giving the final result 

d 2 £dyd 2 b ~2(2vr) 2 



x 



AC f(^) 2 J q dz j d 2 u[l-2S A (u)} 

{ [^((i - ,) 2 + , 2 ) + ™ 2 ] ^FT^F + 4/oV ' ^ } } ' (76) 



Alternatively, we could have used (68) to write 

^ [ d 2 u [ d 2 u'5(u - v!) [1 - 2S A ( 
«r) 4 7 7 



d 2 £dyd 2 bdz 2{2i 

x i / dV2rf 2 ^ 2 e _ ^^ 12 -^2 ) [$ii(a;i,a;2;a;;,4;z) + $22(«i,«2;a;;,a; 2 ;z)] 



+ 2Re J d 2 x 12 d 2 x' 12 e- t£ < Xl2 - x '^^ 12 (x 1 ,x 2 ;x' 1 ,x' 2 ;z)Y (77) 



and plugg in formulas (72)-(74). 

The final result, which is the main result of this paper is 



da q v ^ q qX d<j q ^ X dall^ 

-J27, J27„ ,7„. J2Z, ,7 ' \'°) 



d 2 k\d 2 k 2 dy d 2 bdz d 2 k\d 2 k 2 dy d 2 b dz d 2 k\d 2 k 2 dy d 2 b dz 

where the two terms on the right-hand-side are given by (67) and (75). 

Note that the separation between the contributions o\ and a\\ in (78) is gauge dependent. 
Only the sum of the two cross section has the physical meaning as the heavy quark and antiquark 
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production cross section. In the limit u <C 1/m, it follows that fi — > 1, /o — > 0. If interactions 
are absent, i.e. Sij = 1 this implies that an — > 0, so that o\ (given by (64)) acquires independent 
physical meaning (see the previous section). 

It can be explicitly shown that if we put Sij = 1 for all i,j the cross section vanishes as required. 
Indeed, in this case integrating (67) over x\2 and x' 12 and integrating (75) over u and u' we arrive 
at 



d 2 £d 2 qdyd 2 bdz 
1 /a-\2 



Cf- 



da1^ X 
d 2 £d 2 qdyd 2 bdz 
z 2 (l-z) 2 



.[{{l- z ) 2 + z 2 )£ 2 + m 2 ]+A 



(79) 



(2tt) 2 V 7T / * (z(l - z)q 2 +£ 2 + m 2 ) 2 \ (t 2 + m 2 ) 2 

In order that the g-integrated cross section vanish when Sij = 1, integral over q in (67) must be 
cutoff by q < Q, while integral over u in (76) by u > 1/Q with the same Q. 

The role of the valence quark interactions is illustrated in Fig. 3. It shows numerical calculation 
of charm production with and without valence quark contribution. One observes that collinear 
approximation fails at higher q and lower t. 

10" 4 [ 




FIG. 3: Double inclusive cross section (78) for m = 1.29 GcV (charm) and Q s = 1 GeV. Solid lines include 
contribution of valence quarks, while dashed lines correspond to the collinear factorization of proton. 



V. CONCLUSIONS 

We investigated the collinear factorization of the proton in the process of heavy quark production 
in high energy pA collisions. The factorization holds in the regime Q s <C m, which also implies 
that the gluon saturation effects on the heavy nucleus side are small. However, if Q s > m this 
approximation breaks down and one has to take into account interaction of valence quarks with 
the nucleus. Generally, the collinear factorization of proton is a reasonable approximation only 
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at high relative transverse momentum i of quark and antiquark (high invariant masses M), and 
small total transverse momentum q of the pair (gluon is collinear with the proton), i.e. for hard 
scattering. These observations have important phenomenological implications that will be discussed 
in a separate work. 
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